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Abstract 

We construct bosonic and fermionic matrix-vector models which describe orb- 
ifolded string worldsheets at a limit in which the dimension of the vector space 
and the matrix order are taken to infinity. We evaluate tree-level one-loop or 
multiloop amplitudes of these string worldsheets by means of Schwingcr Dyson 
equations and derive their expressions at the multicritical points. Some of these 
amplitudes resemble or are closely related to those of ordinary multicritical Her- 
mitian matrix models by a constant factor, whereas some differ significantly. 
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1 Introduction 



Large- N matrix models provide us with valuable insights into non-perturbative 
behavior of low-dimensional bosonic strings. (See Ref. jTj and the references 
therein.) This is rendered possible by the observation that the dual of Feyn- 
man diagrams of these models may be regarded as discretised oriented string 
worldshccts and by the tractability of these models at the double scaling limit. 
Recent work has revealed that these models are well suited to the study of 
D-brane dynamics [2121, too. 

There are other important string models besides oriented string theory. For 
instance, one may construct type I superstring theory by an orientifold projec- 
tion of type IIB theory |SJ EI ■ The worldsheets involved are orbifolded and 
respect a Z2 symmetry which interchanges left- and right-movers. Recently, 
we have discovered a family of matrix-vector models which not only serve as 
examples of noncommutative probability of type B [7] but also may be used 
to study models of orbifolded string worldsheets [Sj. The basic ingredients of 
these models are vectors of square matrices of Grassmann numbers. If both the 
vector dimension and the order of the matrices are, loosely speaking, taken to 
infinity, then the Feynman diagrams are the dual of discretised orbifolded string 
worldsheets. It is possible to evaluate the tree-level one-loop amplitudes of the 
simplest of these models. It would certainly be of interest if the calculations 
can be extended to multiloop amplitudes of multicritical matrix-vector models. 
Such calculations are the subject matter of this article. 

Moreover, we will show that there are bosonic counterparts to these fermionic 
models. We will see that the orbifolded string worldsheets that are contructed 
from the bosonic models display some unique characteristics. 

Here is a brief synopsis of this article. In Section|21 we will introduce bosonic 
and fermionic matrix-vector models which describe string worldsheets homco- 
morphic to K 2 /Z2 x Z2. We will derive the tree- level multiloop amplitudes at 
the multicritical points via Schwinger-Dyson equations. In Section |31 we will 
turn our attention to models which describe string worldsheets homeomorphic 
to R 2 /Z2 and use a similar method to evaluate the tree- level one- loop ampli- 
tudes at the multicritical points. Then we will summarise our results and point 
out future directions of this work in Section 0] 

2 Multicritical models of R 2 /Z 2 x Z 2 

Consider a fermionic matrix- vector model whose building blocks are Grassmann 
matrices \& M and "I^ of order N m , where \x may take any integer value between 
1 and N v inclusive and is called a vector index. The action of the model takes 
the form 
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oo N v 
oo JVt, 

3« 



-^E£ E •■•***.-!*/*»)] , (i) 

where c„ and <7„ are constant complex numbers for n = 1, 2, 3, . . . , and so on. 
Like the models we studied in Ref. [Hj, the dual of the Feynman diagrams of 
this model in the double large- N limit in which we take N v to infinity first and 
N m to infinity afterwards may be identified as quadrangulated surfaces of the 
orbifold R 2 /Z2 x Z2. Note that the expression 

Ml,M2 = l 

in Eq. Q may be represented as a pair of Feynman propagators. The term 

N v 



Ni(gi - 1) 



£ [Tr(v& Ml ^ 2 )]- 

Ml,^2 = l 



is put into S 1 / for future convenience. 
Let 

Z f (N m ,N v ) := / c?*id§id* 2 rf*2 • • ■ d^> Nqj d^ Nqj expS/ 



(2) 



be the partition function of this model. The quantities which are of interest to 
us are the connected Green function 

Gf(p 1 ,p2,...,Pn]ki,k 2 ,...,k n ) := lim lim 7V™ +2 "" 2 

iV m — >oo N v -^oo 
N v N v N v 

E E ••• E 

"1,1 ,"2,1, ■■■,U2p 1 , 1 = 1 Kl,2,I'2,3)...,K2p 2 ,2=l fl ,fi ,f2, B , ■ ■ ■ |l^p s , ft = 1 

jV„ JV„ N v 

E E ••• E 

Ml,l ,M2,1 ■■••iM2fc 1 ,1— 1 Ml, 2 ,M2,2 ■•••iM2fc 2 ,2 — 1 Ml,™ I M2,n , • • • ,M2fc„ ,n — 1 



JT ^ ( ''''J' *''"! ' ' ' ''V< 

n 

n [ Tr • • • ^wi,*«w, 



(3) 



conn, 5/ 



where n is any non-negative integer, n is any positive integer, pi, ■ ■ ■ , Pn, fci, 
fc2j ■ ■ ■ j and fc n are also any positive integers, and the subscripts "conn" and Sf 
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v l,l 

(b) 



Figure 1: (a) A Feynman diagram of G s (3;l,2), where s = f or b. This 
diagram corresponds to a vanishing term if s = f . (b) A Feynman diagram 
of G s (2, 1; 1, 2). In both diagrams, the indices mj or Vij , where i and j are 
positive integers, are vector indices. 



4 



tell us that this Green function is connected and that the expectation value is 
evaluated with respect to the action S / , respectively. Terms of some examples 
of Green functions are depicted in Fig. ^ 

There is a bosonic counterpart to the fermionic model. Let Mi, M%, . . . , 
and Mtv„ be complex matrices of order N m . Out of these matrices may be 
constructed a bosonic matrix-vector model whose action is 

oo JV„ 



n=l /ii,^2,— ,M2n = l 

oo JV„ 

-^E^ E [Tr^M^M^M^-.-Aft 



n— 2 Mi,M2,---,A i 2n. — 1 



1 2 



Unlike the fermionic model, those Feynman diagrams of this bosonic model in 
which there is no vertex representing a term whose coefficient is where i is 
any positive integer, do not vanish. Such non-zero Feynman diagrams are also 
invariant under parity transformation. Let 

Z b (N m , N v ) := J dM\dM x dMldM 2 ■ ■ ■ dM^dM^ exp S b (4) 

be the partition function of this model. The physical quantities we would like 
to evaluate are the connected Green functions 

Gb(Pl,P2, ■ ■ ■ ,Pn, ki,&2, • • • , kn) 

defined as in Eq. ® with ^ and the subscript Sf replaced with M\ M, and 
the subscript Sb, respectively. 

2.1 Schwinger— Dyson equations 

We may evaluate the multiloop amplitudes of these matrix-vector models by 
means of Schwinger-Dyson equations. The results are intimately related to the 
ordinary Hermitian matrix model whose action is 

S H := -N m TrV($), 

where $ is a Hermitian matrix of order N m and 

oo 

n*) - E 2n - 



Let 



kn):= lim -^(Tr $ 2 ") Sff (5) 
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be the expectation value of Tr<£> 2n . Consider the trivial equations 

1 



lim lim , 

Zf(N m ,N v ) 



d Nv 

/ V s fill ill * V[r Vl/ . . . * 

QljQ2 ,... iQ2 „_ 1= i 

ex P 5/} = 



13 



and 



lim lim 



w m ^oo n v '00 NlVK>Z b {N m ,N v ) 

N m N v 



E E / dMldM 1 dM%dM 2 ■ ■ ■ dM^dM^ 

A 1 ^. _ 1 J 



i,j=l cc =l 

d 



^^a ij ai,a2,...,Q2n-l=l 

expSf,} = 0, 

where n is an arbitrary positive integer, for the fermionic and bosonic models, 
respectively. They yield the Schwinger-Dyson equation (see Ref. [S] for some 
intermediate steps), 

71-2 



2 H k )Gs(n - 1 - k) + 25 sb j>(n - 1) 

fc=0 

00 00 
-^CkHn + k - 1) -J^gkGsin + k - 1) = 0, (6) 

k=l k=l 

where s = f or b, and S s b (or S s f) is a Kronecker delta function. The first sum 
vanishes if n — 1. Define 

71= 1 S 

and 

71= I S 

as the spectral density functions of the matrix- vector models and of the Hermi- 
tian matrix model, respectively. Then Eq. JfjJ leads to 

^(0= -^(c)-E c ^ 2fe ~ 2 cc) 

I ^ fc=l 
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+ E c 2fe_1 E [*(Oc*+i+i + G(Ofl*+i+i] / 

It is well known (see, e.g., Ref. 9 ) that 

1 r 



E^c 



k=l 



1 1 - - c io 



CO 



where is determined by the integral relation 



W(13) :-- 



qV'(q)dq 
in Jo y/q 2 -[3 



and 



Hence 



M(C,[3) :=£C 2J - 2 E 



fc=0 



(2fc)! 
4 fe (fc!) ; 



Ps(0 



(7) 



(8) 



(9) 



where P S (C) is a polynomial. 

P s (() may be determined by the holomorphic properties of uj s (0- Multiply- 
ing both sides of Eq. © by ±iy/(y/@ - Q{\/P + Q, we obtain 

±^ S (C ± iO)i^/(v^-C)(v^+0 = 



4 + £ CkC ) V(v^-0(v^ + C) + P,(C). 

Thus we get a discontinuity equation 

MC + w) + w s (C - *o)] V(v^-C)(v^ + - 



+ E c ^ 2fe_1 ] *V(v^ - 0(Vp + c). 



fc=i 



As a result, 



.(0 = 



I 



E^ 2fc_1 



fc=i 



Q«(0 



(10) 



where <3s(C) is a polynomial. Since lim^oo w s (C) = 0, Q s {() is constant. Then 
lim^oo £o> s (£) = 1 implies that 



,(C) = i- 
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Evaluating the integral in Eq. I|1(J|) and comparing the result with Eq. © then 
imply 



oo oo 



p s (0 = 28 sb - y ^c 2n + 2 y c 2n Y cfc +™^ 



n=l 



n=0 fc=l 



(2fc- 2)! 
4 fe /c!(fc- 1)!' 



(11) 



We may use Eqs. |J5J| and to expand io s {Q as a power series in l/£ and 
obtain all connected Green functions of the form G s (p). 

2.2 Multiloop correlators 

To obtain other connected Green functions, we apply the formula 

G s (pi,P2, ■ ■ ■ ,Pn] fel,fc2, •• ., k„) 
d 



-2pfi- — G s (pi,p 2 , . . . ,pn-i;h,k 2 , ■ ■ ■ , K) 

OCr, 



— 2k r 



d 

d 9k n 



G s (pi,p 2 , ■ • . ,p fi ;fci,A;2, ■ ■ ■ ,fcn-i)- 



(12) 



Let 



w s (Ci,C2, ■ ■ .,£k;zi,Z2, ••■,z n ) 



E E 



G s (pi,p 2 , ■ ■ ■ ,Ph]ki,k 2 



i fcn) 



>2pi + l >2p 2 +l ^2p s + l 2fei + l 2fe 2 +l 
Pl,P2,...,Ps = l/ci,/c 2 ,...,fc„ = l ^1 ^2 ""Sfi z l z 2 



2fe„ + l 



be the multi-loop generating function of these connected Green functions. This 
may be paraphrased as [H] 



+ 



k=l 

1 d 



c) 



j=i z /t \ u y3/p 



(zl-f3)i \W> ((3)8/3 



9:1 



Ps(0 \ 



(13) 



where (gfrj is the partial differentiation operator with respect to g^- with /3 
held fixed, C^J is the partial differentiation operator with respect to (3 with 

5ii 52i 33j ■ ■ • j and so on held fixed, and W{j3) was defined in Eq. Q. According 
to Ref. 0, 



f> JLfJL 



3 = 1 Z k 



dg 3 J p (zi-P)i \W'(0)d/3 /;i 
1 d 1 " +1 /i(/3) p 



_J_d_ 

W'((3) dp 



W'{(3) 



W'{p) dp 



w>(j3) (4-p)\ 
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if n is a non-negative integer and h{(3) is a function which depends only on (3 
but not gi, g 2 , 53, ■ ■ ■ , and so on. As a result, 



n-l 



1 



(9 



W(/3) 9/3 



fe=l 



(2fc-l)! 
4 fe fe!(fc-l) 



— 5* 
1)1 V 







for any positive value of n. 

In addition, Eq. i|12|) implies 



W s (Cl,C2) 



2C1C2 



2Ci 2 Cl 



/3(Ci 2 + CD 



(C1 2 - CI) 2 (C 2 - C 2 2 ) 2 v / CT^n/CT^' 



(14) 



Note that u> s (Ci, C2) is independent of whether the model is bosonic or fermionic 
and is independent of c±, c 2 , C3, . . . , and so on. Thus we conclude from Eqs. I|14(l 
and JTJJl that 

W s (Cl,C2, ■ ■ -,Cn,Zl,Z2, ■ ■ -,Z n ) = 



if n > 3. In other words, 

G s (pi,p 2 , ■ ■ ■ ,f>«; fci, &2, 



, fc„) = 



if n > 3. In terms of string worldsheet, this means that there can be only 
two boundaries which are invariant under parity transformation. Note also that 
Eq. (|14|l differs from the two-loop correlator of any complex matrix model by 
a factor of 4 only Since w s (Ci,C 2 ) depends on g\, g 2 , 33, . . . , and so on 
indirectly via /3 only, we could apply a formula similar to Eq. Ijl3|l to obtain 
other generating functions: 



■ j Zn) — 



" 1 d ' 


n-l 


1 (i (3 


_W'(f3) df3_ 




20W'(P) (Ci 2 -/?)* (C 2 2 -/5)3 



n 







(15) 



for any positive value of n. These multiloop generating functions differ from 
those of complex matrix models merely by constant factors of 2 ,l+2 [5]. They 
are basically symmetry factors of the Feynman diagrams. 

2.3 Multicritical point 

Following Ref. 10j . we approach the m-th multicritical point by fine-tuning 
the coupling constants in such a way that there exists a real number z c which 
satisfies 

W(0*) = W (k) ([3*) = 
for k = 1, 2, . . . , and m — 1, and 
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Then 

where 7 is a complex constant, for /3 close to z*. Let 

(f = /3* + aipi, z\ = f3* + air l , and (3 = /3* - aVX 

for any positive integer i, where a is the cut-off length, A is the renormalised bulk 
cosmological constant, and 7Tj and ipi are renormalised boundary cosmological 
constants for any value of i. Then 



1 d 



d 



W>(/3) dp m 7 (-a)™A^ OA 
and we may conclude that the renormalised tree-level one-loop amplitude is 



/aw s (V>i) 



OA^ + V 00 /?"r 2 V" (2fc-2)! I 



and the renormalised tree-level multi-loop amplitudes are 

a ("H-i)n+i a , i (^ i;7ri)7r2j ... >7rn ;A) = 

/i^mn+ti+l9Ti-2on i _ -sr^ 00 (2fe-l)! ofe 

I JJ z P* 1 Z^fc=i iJ^PIfc^T)!^* 



<9 



n-l 



1 



1 



n 



i 



for n > 1, 

a 2 u; s (V>i,^2) 
and 



A 2 ?-^ 9A; + VA)i (TTi + VA)i_ 

ty'l + V^2 + 2y/A)0. 



2/3* 



+ 



(V-l-^) 2 (V'l-^) 2 (^1+^)5(^2+^)^ 



a (m+i) " +2 ^(^i,^ 2 ;^^2,...,7r„;A) 
(-l)" m +"/3™ +1 / 1 9 



4to™7™ 



<9A 



n 



_A^(Vi + VA)i(^ 2 + VA)i + VA)i 



for n > 1. 
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3 Multicritical models of R 2 /Z 2 

Let us turn our attention to multicritical models of the quantum orbifold R 2 /Z2. 
The action of the bosonic version is 

oo N v 

-N^cn E Tr (m^M^M^M^.-.M^M^ 

n=l Hi, 112, ...,/i2n=l 

■M^Ml^M^Ml^ ■ ■ ■ M^Ml 

oo N v 

■Tr (M^Ml^M^Ml^ ■ ■ ■ M, 2 Ml) , 
whereas the action of the fermionic version is [5] 

5/ := -N m V^E Tr *^^ 

g Tr^^jTr^^) 

Ml ,M2 = 1 

oo N v 
n=l Hi,Ha,— ,/i2n=l 

'^P2n ^M2n-1 ^^271-2 ^^2^-3 ' ' ' ^A^^Vl) 

oo W„ 

— TV 2 V* — Trfvfrxl/vi/fl/ . . . \J/ * ") 

•Tr (^^-^M^^Wa • • • ■ 

Note that the second terms in these actions may be represented by a pair of 
Feynman propagators. The partition functions of the bosonic and fermionic 
models are defined as in Eqs. (@J and J3J, respectively, with Z replaced with Z 
and 5 with S. For the bosonic model, the connected Green functions which we 
would like to study take the form 



G b (pi,p 2 ,.-.,Pn,ki,k 2 ,.-.,k n ):= lim lim iV, 

N m —>oo N v —>oo 



n+2n-2 



11 



V(,l 




Figure 2: A Feynman diagram of Sb(3, 1, 1; 2). 



E E - E 

"1,1 ,"3,1 I'-'.^P!, 1 = 1 ^l,2,^2,2,---,"2p 2 ,2 = l )V2,ft i • •M 1/ 2p a ,&=l 

AT„ AT„ JV„ 

E E ••• E 

Ml, 1 1^2,1 ,...,(12*! ,1=1 Ml, 2 ,M2,2,---,M2fc 2 ,2 = 1 /*l,n ,M2 ,n , ■ ■ ■ ,M2fc„ ,n=l 

( TT Tr f AfJ M„, . M I M Vi ■■■ Ml M Vo . 

\i=i 

n 

nTr f M,t M u , .M,t M u . . • • • M,t M„,, 
V Mi,i M2,s )i3,, M4,> M2fc,-i,. f*afc<,« y 

»=1 

■^(M^Mt ^^Af^^Ml 2fci _ 3 , i ---M M2ji Mt ii )\ , (16) 

V 1 1 ' / / conn, St 

where n is any non-negative integer, n is any positive integer, and p%, P2, ■ ■ ■ , 
Pn, ki, k2, . . . , and k n are any positive integers; for the fermionic model, the 
connected Green functions which we would like to study also take the form 
in Eq. (|16l) with M, M', and Sb replaced with 'J, 'J, and Sf, respectively. A 
Feynman diagram representing a term in a connected Green function is depicted 
in Fig. El 
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3.1 Schwinger-Dyson equations 

To evaluate the connected Green functions at the double large- N limit, let us 
start with the trivial equation 

lim lim 



N m N v 



J dM\dM l dMldM 2 ■ ■ ■ dM f N JM Nv 



d 



—{ £ (Ml n _ 1 M a2n _ 2 ---Ml 1 



^M ao ij a 1 ,a 2 ,...,a 2n _ 1 = l 



■ Mi x M a2 ■ ■ ■ Ml 2n i M ao ) exp S b } = (17) 
for the bosonic model or 



1 N m N 

lim lim == / d^idty 1 d^ 2 d^2 ■■■ d^ N d^ N 

N m ^N^ooNlV7r v Z' f {N m ,N v ) .j^^J N * N * 



d Nv 

1 > (\tl \J/ . . . \p \[; \i; ...X]/ \Tf ^ 



u J Qi ,a2 ,. . . ,Ct2n — 1 = 1 

expS*/} = (18) 

for the fermionic model. Both Eqs. I|17l) and l|18[l lead to the Schwinger-Dyson 
equation 



Y4>(n~ k)G{k - 1) 

k=l 

oo k oc 

-^c fc ^G(Z-l)G(n + fc-0-E^( n + fc-1 ) = O ' ( 19 ) 

fc=l 1=1 k=l 

where n is any positiver integer, G(n) stands for Gb(n) or Gf(n), and 4>{n) was 
defined in Eq. JSJ . Hence the connected Green functions of the bosonic model 
at the double large- N limit are identical to those of the fermionic model. 
Let 

^ G{n) 

n=l ^ 

be the spectral function of these matrix-vector models. It then follows from 
Eq. (|19l) and the well-known expression for the spectral function (£) of the 
ordinary Hermitian matrix model that 

- (0 = 2Ql(C) - (20) 

Q 2 (c)+M(c,i3)Ve^p 
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where (3 and M((,(3) were denned in Eqs. @ and |[SJ, 



(21) 



Qi(o :=E^ 2 E«^) +E^E 6 (0 E GH^i+m+i, (22) 



fe=l 2=0 fc=0 i=0 m=0 

and 



Q 2 (C) :-£. 9fc C 2 ' £ - 1 +2£c 2fc+1 E (5 ( / ) c ^+i- (23) 

fe=l k=0 1=0 

As usual, we assert that the values of the connected Green functions in Eqs. 12 I'll 
and (|23[1 are determined by the requirement that uj be holomorphic on the whole 
complex plane except the branch cut — y 7 /? < 3?(C) < V/3 an d 3(C) = 0- 

3.2 Some multicritical points 

A convenient choice of the m-th multicritical point is to select a non-zero value 
of g m , adjust the values of g±, g2, ■ ■ ■ , and g m -i such that 

M(C) = g m (C 2 -P) m -\ 

and adjust the values of ci, C2, . . . , and c m such that 

02(C)=3 m C(C 2 -/3) m - 1 . (24) 
Moreover, g n = Cn = if n > m > 2. It then follows from Eq. (121 III that 



2 Qi(C)(C-v / C T ^) 

5m /3(C 2 - /3) r 



-(0 = „ .;; 2 Vm-i • (25) 



The holomorphic property of U)(Q then dictates that the zeros of Qi(£) coincide 
with the zeros of the denominator on the right side of Eq. Q25J1. As a result, at 
the m-th multicritical point, 

Qa(C) = A(C 2 - /9) m " x ; (26) 
the constant A may be determined by the condition that 

lim C<y(C) = 1. 

This yields A = g m . As a result, 



(0 = |(c-Vc^) 
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at the m-th multicritical point. 

A convenient way to approach the m-th multicritical point is to keep the 
ratios gi : gj , gi : Cj , and Cj : Cj , where i and j are positive integers less than or 
equal to m, fixed. Then in Eq. (|2(J[1 . only £ and (3 deviates from their critical 
values v 7 ^* an d P*, respectively, whereas Ck, and G(fc), where k is any 
positive integer not larger than m, are fixed. Let 

C 2 = P* + an and /3 = /3* - aVX, 

where a is the cut-off length, and tt and A are the boundary and bulk cosmolog- 
ical constants, respectively. Then Qi((), Q2((), and M(£) are of order a m ~ l , 
whereas \/ ( 2 — f3 is of order yfa. Recall that M(£)a/£ 2 — f3 is, up to a pro- 
portionality constant, also the singular part of the spectral function of ordinary 
Hermitian matrix models. Hence we conclude from Eqs. (|21J) . (|24(l . and l|26|) 
that we may multiply u) by ^fa to obtain the renormalised tree-level one-loop 
amplitude which, up to a constant factor, is 

^ / renormalised 1-loop amplitude of \ 

=- an ordinary Hermitian matrix model 

\ at the m-th multicritical point / 



4 Conclusion and Outlook 

We may study quantum orbifold geometry by means of bosonic or fermionic 
matrix- vector models. As for the quantum orbifold K 2 /Z2 x Z2, the bosonic 
model differs from the fermionic model in the sense that Feynman diagrams 
with no Ci-vertices, where i is any positive integer, contribute to the Green 
functions of the bosonic model only; they have no contribution to those of the 
fermionic model. If in an orbifolded worldsheet there is only one boundary 
which is invariant under parity transformation, then its multiloop amplitude 
is significantly different from that of an ordinary worldsheet. Nonetheless, if 
there are two boundaries which are invariant under parity transformation, then 
its multiloop amplitude is the same as that of an ordinary worldsheet up to a 
symmetry factor. 

As for the quantum orbifold M 2 /Z2, the bosonic and fermionic models are 
equivalent to each other at the double large- N limit. The renormalised tree- 
level one-loop amplitude at an m-th multicritical point differs from that of an 
ordinary Hermitian matrix model by a factor inversely proportional to 7r m_1 . 
Nevertheless, it may be possible to identify other m-th multicritical points at 
which the quantum orbifold may behave differently It would also be of interest 
to obtain more explicity expressions for higher loop amplitudes of this quantum 
orbifold. Furthermore, exploring the double-scaling limit of these matrix-vector 
models would give us valuable information on the non-perturbative behavior of 
unoriented string theory. 
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